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We report on the synthesis, scanning tunneling microscopy (STM) and theoretical modehng of 
the electrostatic and transport properties of one-dimensional organic heterostructures consisting of 
contiguous lines of CF3- and OCHa-substituted styrene molecules on silicon. The electrostatic fields 
emanating from these polar molecules are found, under appropriate conditions, to strongly influence 
electrical conduction through nearby molecules and the underlying substrate. For suitable align- 
ment of the OCH3 groups of the OCHa-styrene molecules in the molecular chain, their combined 
electric fields are shown by ab initio density functional calculations to give rise to potential profiles 
along the OCHs-styrene chain that result in strongly enhanced conduction through OCHs-styrene 
molecules near the heterojunction for moderately low negative substrate bias, as is observed exper- 
imentally. Under similar bias conditions, dipoles associated with the CF3 groups are found in both 
experiment and in theory to depress transport in the underlying silicon. Under positive substrate 
bias, simulations suggest that the differing structural and electrostatic properties of the CF3-styrene 
molecules may lead to a more sharply localized conduction enhancement near the heterojunction at 
low temperatures. Thus choice of substituents, their attachment site on the host styrene molecules 
on silicon and the orientations of the molecular dipole and higher multipole moments provide a 
means of differentially tuning transport on the molecular scale. 

PACS numbers: 31.70.-f, 68.37.Ef, 68.43.-h, 73.63.-b 



I. INTRODUCTION 

During the last decade a great deal of research has 
focussed on electrical conduction through individual 
molecule d^ l ^ l '^l Many molecules are electrically polarized 
due to chemical charge transfer between unlike atoms. 
This results in electric fields that should infiuence elec- 
trical conduction through such molecules and through 
molecules in their vicinity. However, there have been 
few direct experimental investigations of such effects in 
the context of molecular-scale nanoelectronics. Recently 
experimental studies of the effects of charged chemical 
species attached to a molecule on electrical conduction 
through the same molecule have been reportecP^J. The 
presence of a charged dangling bond on a silicon surface 
has been observed to affect electrical conduction through 
nearby moleculed^and, conversely, transport through ad- 
jacent silicon atoms has been found to be perturbed by 
dipole fields due to molecules located elsewhere in the 
Si 7x7 celP. However, no experimental work directly 
probing the effects of electric fields emanating from po- 
lar molecules on electrical conduction through other indi- 
vidual molecules has been reported to date. This topic is 
explored experimentally and theoretically in the present 
article The influence of these electric fields on elec- 
trical conduction through the underlying semiconduc- 
tor substrate to which the molecules are bound is also 
examined .^^ 

The specific systems that we study here are one- 



dimensional (ID) organic heterostructures consisting 
of chains of substituted styrene molecules grown on 
hydrogen-terminated Si substrates by the self-directed 
growth mechanism described by Lopinski and coworker^ 
originally for styrene on H:Si(100) 2x1 but subsequently 
applied to a wide range of molecules. The substituents in 
the present work are CF3 and OCH3 groups, one of which 
replaces a single hydrogen atom bound to the aromatic 
carbon ring of each styrene molecule. In each case the 
substituent is in the "para" position. Each heterostruc- 
ture consists of a line of surface bound OCHs-styrene 
molecules joined end-to-end (at the heterojunction) to 
a line of CFa-styrene molecules. The OCH3 group do- 
nates electrons to the aromatic ring of the styrene and 
is therefore positively charged while the CF3 group with- 
draws electrons and is negatively charged. However, elec- 
tric polarization also occurs within the CF3 and OCH3 
groups themselves and we find the resulting molecular 
multipole fields to also play an important role in electri- 
cal conduction in these systems. 

In the experimental work presented here electrical con- 
duction between the silicon substrate and a tungsten 
scanning tunneling microscope (STM) tip via the molec- 
ular heterostructures was measured at room tempera- 
ture in ultra-high vacuum under a variety of bias con- 
ditions. Strongly enhanced conduction was observed 
through a group of several 0CH3-styrene molecules near 
the heterojunction at low and moderate negative sub- 
strate bias. Under positive substrate bias similar pro- 
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nounced enhancement of the molecular conduction near 
CFs-styrene/OCHs-styrene heterojunctions was absent. 
In related heterostructures, CFs-styrene lines in a side- 
by-side configuration were found to locally depress trans- 
port originating from filled states in the underlying sil- 
icon. As will be explained below these qualitative dif- 
ferences can be understood as arising from the differ- 
ent structural, electronic and electrostatic properties of 
OCHa-styrene and CFs-styrene molecules on silicon. 

The OCHa-styrene/CFa-styrene heterowires were 
modeled theoretically by determining their ground state 
electrostatic potential profiles (at zero applied bias) by 
means of ab initio density functional theory-based com- 
putations, and then employing semi-empirical tight bind- 
ing models together with the ab initio electrostatic pro- 
files, solution of the Lippmann-Schwinger equation and 
Landauer theory to calculate the electric current between 
the tungsten STM tip and the silicon substrate under bias 
via the molecular heterostructure. We note that while 
the use of density functional theory for performing trans- 
port calculations in molecular wires lacks a fundamental 
justificati on and is increasingly being que stioned in the 
literature,I^Eo|i i | i^ l " | i4 | ^^ | i^ | i^ l ^« | i" l ^" i ^i | ^" lthe use of den- 
sity functional theory to calculate ground state electro- 
static potentials is justified (at least in principle) by a 
lemma proved originally by Hohenberg and Kohn.'^'^ The 
usefulness of such density functional theory-based elec- 
trostatic potential calculations for the present system is 
supported also by the fact that the present theory is able 
to explain the experimentally observed phenomena out- 
lined in the preceding paragraph. 

The theoretical work that is presented here demon- 
strates that electric fields due to intramolecular 
charge transfer within OCHa-styrene and CFa-styrene 
molecules can (for appropriate molecular geometries) re- 
sult in enhanced electrical conduction through certain 
molecules in chains of OCHa-styrene and CFa-styrene 
molecules on silicon, consistent with the experimental 
STM data: These electric fields shift the energies of the 
molecular HOMO and LUMO states of like molecules in 
different parts of the chain by differing amounts with 
the result that resonant conduction begins at lower bias 
voltages on some molecules in the chain than on oth- 
ers. The differences in the observed behavior under pos- 
itive and negative applied bias and between the OCH3- 
styrene and CFa-styrene molecules are accounted for in 
terms of this electrostatically modulated resonant trans- 
port mechanism. 

The present theory predicts that conduction through 
individual molecules in such systems can be changed by 
orders of magnitude by varying the conformations of (and 
hence the arrangement of charges in) molecules in their 
vicinity. This raises the possibility of molecular switches 
of a new type that depend for their operation not on con- 
formational changes in the molecular wi re carry ing the 
current (as has been discussed previousl> *^^ * ^^*^^ *) but on 
conformational changes in other nearby molecules that 
result in changes in the energy level structure of the 



molecular wire. An advantage of a conformational switch 
based on this different principle is that the molecular wire 
that carries the current need not be a moving part in the 
device. The shorter range of electrostatic fields due to 
molecular dipoles (and higher multipoles) than those of 
electric monopoles such as the charged dangling bond 
studied in Ref. ^ is also advantageous, making it possi- 
ble in principle to achieve higher densities of molecular 
switches controlled by molecular dipoles in potential de- 
vice applications. The greater chemical stability of a po- 
lar molecule than a charged radical is also an important 
advantage. Furthermore manipulating the orientations of 
molecular dipoles may be an attractive alternative^^' to 
inducing switching by charging and discharging atomic 
or molecular-scale constituents of the ultimate nanoscale 
electronic devices. 

This article is organized as follows: The experimental 
methodology and results are presented in Section [TTj The 
theoretical model is explained and justified in Section [ITT] 
A discussion of the relevant energy level ordering in these 
systems is presented in Section [ill B[ The different struc- 
tures that the molecules may assume when bound to sili- 
con are described in Section UlI CI Our theoretical results 
for a particular conformation of a OCHa-styrene/CFa- 
styrene molecular chain on silicon are presented in Sec- 
tion |TV] for positive and negative substrate bias, together 
with some comments as to how these theoretical results 
may relate to the experimental data. In Section |V] we 
clarify the relationship between structure, electric po- 
tentials and transport in the OCHa-styrene/CFa-styrene 
heterostructures by considering systematically other ex- 
amples of possible molecular chain geometries. Theoret- 
ical results demonstrating that the current enhancement 
near the OCHa-styrene/CFa-styrene junction is specif- 



ically an electrostatic effect are reported in Section VI 



Simulations of heterostructures that include single and 
triple rows of CFa-styrene molecules and of the influence 
of these molecules on the electrostatic potentials in the 
underlying silicon and on electron transport are reported 
in Section |VII[ Further discussion of the relationship be- 
tween the theory and experiment is presented in Section 

Ivml 



II. EXPERIMENT 

STM experiments were performed under vacuum on 
hydrogen-terminated Si (100) 2x1 surfaces. Samples 
were cleaved from arsenic-doped (resistivity < 0.005 
ricm) Si (100) oriented wafers, and mounted into molyb- 
denum sample holders. Samples were loadlocked into a 
vacuum system (background pressure < lxlO~^° Torr) 
and degassed at 700 °C for 8 hours. The samples were 
flash annealed to 1250 °C to remove the surface oxide and 
re-order the crystal surfaces. During annealing, heating 
was temporarily suspended if the system pressure ex- 
ceeded 4x10"^" Torr. Clean crystalline surfaces were 
routinely produced with defect densities below 5%. 
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After cool down and inspection in the STM, the silicon 
crystals were transferred to the preparation chamber for 
hydrogen termination. Molecular hydrogen was leaked 
into the system (1x10^^ Torr) and a hot tungsten fila- 
ment (~1600°C) positioned 10 cm from silicon sample 
(T — 300 °C) dissociated the molecular gas into reac- 
tive atomic hydrogen. A 13 minute exposure produced 
samples with a quasi-saturated 2x1 silicon monohydride 
surface .'^Sl 

One dimensional molecular organic heterowires were 
grown under vacuum using the self-directed growth 
mechanism reported by Lopinski and coworker^ for 
styrene on H:Si(100) 2x1. For styrene (consisting of 
an aromatic carbon ring bound to a vinyl group), self- 
directed growth results from a chain reaction between the 
vinyl group of the styrene and a surface dangling bond 
(i.e. silicon radical exposed by a missing hydrogen atom 
on the surface). In forming a chemical bond with the 
exposed silicon atom, the terminal carbon atom on the 
vinyl group breaks one of its double bonds to the adjacent 
carbon atom, leaving the unsatisfied bonding electron on 
the adjacent carbon atom to abstract a hydrogen atom 
from an adjacent Si dimer. The newly formed dangling 
bond is then free to repeat the process by reacting with 
another styrene molecule. Multiple reactions lead to well 
ordered one dimensional chains of styrene bound along a 
given side of a silicon dimer row.lSl 

This self-directed growth mode has been observed 
to occur for various alkene and carbonyl contain- 
ing molecules including functionalized styrenePSEoMEIi 
making the formation of heterowires a simple matter 
of sequentially dosing the H: silicon surface with the de- 
sired chemical precursors.-^' The fortuitous alignment of 
the aromatic rings in these structures make them ide- 
ally suited for probing transport effects resulting from 
overlapping tt and tt* states along the molecular chain 
axis. The introduction of a chemical discontinuity at the 
heteroj unction provides a means of isolating and study- 
ing effects resulting from specific intermolecular interac- 
tions. Electron donating and withdrawing substituent 
groups are of particular interest as they modify the en- 
ergy alignment and spatial distribution of tt and vr* 
states in host aromatic molecules. These effects lead to 
systemati c var iations in reaction rate^^ and ionisation 
potentiald^HH in substituted aromatic compounds. Cor- 
responding effects on single molecule transport are under 
current investigation.- 

Substituent effects in STM imaging of 4- 
methylstyrene/styrene heterowires were reported in 
Reference ^El While the substituted methyl group had 
a strong (bias dependent) influence on the differential 
molecular height and corrugation resolved on either 
side of the heterojunction, the imaging characteristics 
of the heterowires could be understood qualitatively 
without considering the effects of intermolecular elec- 
trostatic interactions on the electronic structure of 
the heterowires; features reflecting such interactions 
were not observed in the experimental data. The 



present work revisits the earlier heterowire line-growth 
experiments with more strongly perturbing substituent 
groups P^l^ 4-methoxystyrene (OCHs-styrene) , and 
4-trifluoromethylstyrene (CFs-styrene) . 

H-termination of the Si(lOO) samples, growth of the 
molecular heterowires and their imaging with an Omi- 
cron STMl were all carried out within the same vac- 
uum system. Dissolved atmospheric gases were removed 
from the substituted styrene precursors using multiple 
freeze-pump-thaw cycles prior to introduction (via leak 
valve) into the vacuum system. Before and after line 
growth, samples were imaged in STM using electrochem- 
ically etched tungsten tips, cleaned by electron bombard- 
ment, and fleld emission. Bias dependent STM imag- 
ing of the structures was performed in constant-current 
mode (tunnel current flxed at 40 pA). Surface coordi- 
nates belonging to STM images in Figures [l] to |3] were 
rescaled after acquisition to yield an orthogonal inter- 
row (0.768 nm) to inter-dimer (0.384 nm) periodicity of 
2:1. Feature heights in the constant-current STM images 
were determined in relation to observed terrace height 
separation on the Si(lOO) surface (0.136 nm). Images 
of CFa-styrene/OCHs-styrene heterowires were acquired 
on multiple surfaces and studied using several different 
tips in two separate Omicron STMl systems. 

Figure [l] shows the growth and bias dependent STM 
imaging of two CF3-styrene/OCH3-styrene heterowires 
on an H:silicon sample. After imaging the unreacted 
H:silicon surface (not shown), the STM tip was retracted 
^1 micron from the surface. CFs^styrene was introduced 
into the chamber at a pressure of 1x10^^ Torr for ^^10 
seconds. Sample exposure was nominally lOL (IL — 10^^ 
Torr sec). The H-terminated silicon dimer rows run di- 
agonally between the upper left and lower right hand 
corners of the image frames. The elongated white (ele- 
vated) features running along the dimer rows correspond 
to regions where molecules have reacted with the surface. 
The black (depressed) features appearing most notably 
in the centers of Figs [T] (a) and (b) result from missing 
atoms (Si vacancies) in the silicon surface. 

Fig|lja) shows a 26nmx26nm region of the sample fol- 
lowing the lOL exposure of CFa-styrene. Sample bias 
(Vs) was —3.0V. White arrows label the reactive dan- 
gling bonds at the end of two CFa-styrene hne segments. 
Due to slight tip asymmetry, CFs-styrene bound to ei- 
ther side of their host dimers image with slightly different 
corrugation. Comparison with images of the unreacted 
H:silicon surface (not shown) show the upper and lower 
CFa-styrene line segments are chemically bound to the 
right-hand and left-hand sides of their respective dimer 
rows. 

Figjljb) shows the same region of the crystal {Vs = 
—3.0V) following a lOL exposure of OCHa-styrene (1 
X 10^^ Torr exposure for ^10 seconds). The termi- 
nal dangling bonds belonging to the upper and lower 
CFa-styrene line segments in (a) have reacted with 
OCHs-styrene forming two CFs-styrene/OCHs^styrene 
heterowires ('1' and '2'). At —3.0V, the tip Fermi-level is 



4 



below the highest band of occupied molecular tt states for 
the OCHa-styrene as evidenced by the fact that at this 
bias voltage the experimental STM height profile of the 
OCHs-styrene has saturated as can be seen in Fig. 1(d): 
This saturation indicates that the number of 0CH3~ 
styrene HOMO states contributing to the STM current 
is no longer increasing with increasing bias so that the 
STM tip Fermi level must be below the highest band of 
OCHa-styrene HOMO levels. However, the tip Fermi- 
level remains above the occupied molecular tt states in 
the CFa-styrene at this bias because the molecular CFa- 
styrene HOMO is well below the OCHg-styrene HOMO 
as evidenced by gas phase molecular calculations'^^ and 
by the present theoretical work for these molecules on 
silicon. This is also consistent with the CF3- styrene 
molecules imaging with reduced height (less bright) in 
comparison to the OCHg-styrene in Figjl] 

Fig|l|c) shows the same region of the crystal imaged 
with Vs = -1.8 V. At this bias, the tip Fermi-level is in 
the vicinity the highest occupied molecular orbitals (tt 
states) in the OCHa-styrene line segment; this conclu- 
sion is based on comparisons between these experimental 
data and the results of our theoretical modeling as is dis- 
cussed in Section [IVB[ While the OCHa-styrene contin- 
ues to image above (brighter than) the CFa-styrene, het- 
erowires 1 and 2 image with enhanced height (brighter) 
above the OCHa-styrene molecules situated near the het- 
erointerface. The OCHa-styrene molecules in heterowire 
1 (close-up in inset) near the terminal dangling bond also 
image with increased height. 

Fig|T|d) shows a series of 0.4nm wide topographic 
cross-sections extracted from heterowire 1 along the 
trench running between its attachment dimer row (la- 
belled with red circles in FigQc) inset) and the vacant 
H-terminated dimer row to its right. Also included are 
curves from images (not shown) acquired at intermediate 
sample biases. The topographic height envelope for the 
heterostructure extends between 1 nm and 9.5 nm along 
the abscissa. The height maxima associated with the ter- 
minal dangling bond, and the molecular heteroj unction 
are located at 2.3 nm, and 6.4 nm, respectively. Post- 
poning discussion of the interfacial height enhancement 
for the time being, the bias dependent height response of 
the OCHa-styrene line segment near the terminal dan- 
gling bond is much like that reported in ReflHlfor a single 
homomolecular styrene line segment: At elevated nega- 
tive substrate bias, (|T4| ^ 2.4 V) the highest band of 
occupied tt states in the OCHa-styrene lies above the tip 
Fermi-level and the OCHa-styrene images with roughly 
uniform height. As the bias decreases in magnitude, 
these occupied tt states of the OCHa-styrene approach 
the tip Fermi-level (and eventually begin to drop below 
it), and the molecules image with reduced height above 
the surface. Molecules in the vicinity of the negatively 
charged dangling bond on the degenerately doped n-type 
surface exhibit a spatially dependent reduction in ioni- 
sation potential due to the electrostatic fiel d em anating 
from the negatively charged dangling bond.'^'^ Thus at 



low bias, molecules nearest the terminal dangling bond 
image with increasing height.!^ 

In stark departure from the imaging characteristics of 
homomolecular styrene chains reported in Ref.6, and ab- 
sent from the images of the CHa-styrene/styrene het- 
erowires studied in Ref|35] is the height enhancement 
at the heterojunction resolved at low bias. The bias- 
dependent height response of the OCHa-styrene near the 
heterojunction is similar to that near the terminal dan- 
gling bond just described. At elevated bias, the inter- 
facial OCHa-styrene images with nearly constant height 
along the bulk of the homowire segment. As \Vs\ de- 
creases, however, the height of the interfacial OCHa- 
styrene (4-5 molecules closest to the molecular interface) 
does not decay as rapidly as in the mid-section of the 
OCHa-styrene line segment. At the lowest fiUed-state 
bias studied for this particular heterowire, the interfacial 
OCHa-styrene molecules image ~0.05 nm higher than 
OCHa-styrene molecules situated 5-7 dimers away from 
the heterojunction. The typical noise in the presented 
images is < 0.004 nm. 

Figj2] shows filled and empty-state imaging of a clus- 
ter of 3 heterowires. These heterowires were identified 
on a region of the crystal where imaging of the H:silicon 
surface before line growth was not carried out. The chem- 
ical identities of the line segments were confirmed at the 
end of the imaging sequence by dosing with CFa-styrene, 
and comparing the imaging characteristics of the newly 
reacted line segments (not shown) with the nearby het- 
erowires. Characteristics in these images (in combination 
with those in FigjlJ are representative of the range of tip- 
dependent imaging contrast encountered over the course 
of the experiments. 

Figj2]ja) shows a fiUed-state constant current STM im- 
age of the 3 heterowire cluster with Vg — —2.6V. As in 
Figjljb), the OCHa-styrene line segments image taller 
(brighter) than the CFa-styrene line segments consis- 
tent with their relative ionisation potentials. In Figj2]^b) 
{Vs = —2.0V), however, the OCHa-styrene away from 
the interface images lower (darker) than the CFa-styrene. 
This contrast reversal (tip dependent) is often encoun- 
tered at low bias in fiUed-state imaging. In all instances, 
whether or not this low bias contrast reversal occurs, 
height enhancement of the interfacial OCHa-styrene re- 
mains prominent and corresponds to the height of OCHa- 
styrene far from the interface at greater magnitude fiUed- 
state bias (i.e. the interfacial OCHa-styrene images as 
though accompanied by a localised increase in effective 
tip-sample bias magnitude). 

Figs[2]^c) to (f) show empty-state images for the 3 het- 
erowire cluster. In Figgc) at Vs = +2.6Y, the CFa- 
styrene images above (brighter than) the OCHa-styrene. 
As is discussed in Section IIVBI this result is consistent 
with the increased electron affinity of the CFa-styrene 
(a result of the highly electronegative F atoms). As the 
TT* state for the CFa-styrene is much lower than the cor- 
responding TT* states belonging to the OCHa-styrene, a 
greater number of tip states can tunnel into the CFa- 
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styrene (and therefore, the STM tip must travel further 
away from the surface above the CFa-styrene to maintain 
a constant tunnel current across the heterostructure) . In 
Figj2]^d), Vs = +2.0V and the tip Fermi-level is presumed 
to lie below both the OCHs-styrene (as before) and CF3- 
styrene tt* states. The tunnel current is dominated by 
carrier injection into the silicon conduction band, and the 
molecules image (non-resonantly) with similar contrast. 

Figs |2]^e) and (f) display empty-state imaging data 
acquired at Vs — +2.0 V following spontaneous tip 
structural changes. In this imaging mode, the OCH3- 
styrene images above (brighter than) the CFa-styrene. 
This imaging behaviour when apparent, often extends to 
greater bias (up to Vs — +3.2 V). In addition to this 
empty-state contrast reversal, changes to the molecular 
corrugation within the homomolecular line segments are 
apparent (particularly in terms of the appearance of the 
OCHs-styrene between frames (d) and (f)). It is clear 
in these experiments that the tip density of states plays 
a considerable role in determining the contrast (in terms 
of overall height and/or molecular corrugation) observed 
between the OCHs-styrene and CFa-styrene line seg- 
ments. Tip-dependent imaging contrast in STM of or- 
ganic molecules has been documented previously. ^^'^ 

While the empty-state imaging displays considerable 
contrast variation, of particular significance is the ab- 
sence of notable interfacial structure. An exception to 
this was observed in a small portion of low bias empty- 
state images which revealed height enhancement for a 
single interfacial molecule. The height enhancement was 
between 0.01 to 0.02 nm, and smaller than the 0.02 nm 
to 0.14 nm height enhancement observed for the consid- 
erably broader 5 molecule wide) fiUed-state interfacial 
features resolved in Figs [TJc)-(d), |2jb), and others (not 
shown). While such a narrow empty-state feature ap- 



pears in the simulation work presented in Sections IV 



and |V] additional observations will be required before a 
detailed comparison with theory can be undertaken. 

Fig. [3] shows STM imaging results for a CF3- 
styrene/OCHa-styrene heterostructure with multiple 
CFs-styrene chains in a side by side configuration. Fig. 
[Sf^a) shows a region of the H:silicon surface following a 
10 L exposure of CFa-styrene {Vs = —3.0 V). The ar- 
row points to the reactive dangling bond at the end of 
the longest CFa-styrene line. The ★ marks a short dou- 
ble chain of CFa-styrene that has grown beside the long 
CFs-styrene chain. Figures [3](b) to ^d) show the same 
region following a 10 L exposure of OCHa-styrene. The 
end of the long CFa-styrene chain has been extended 
by approximately 7 molecules of OCHs-styrene. Figures 
[ajb) {Vs = +2.0Y) and [^c) {Vs = -3.0V) image the 
single and triple CFa-styrene segments with comparable 
height. In Fig. |3jd), Vs has been reduced to —2.0 V and 
the region with the triple CFa-styrene lines images below 
(darker than) the single file chain of CFs-styrene. Fig. 
^e) shows topographic cross sections along the CF3- 
styrene/0CH3-styrene heterowire. From Vg = —3 V to 
Vs = —2 V, the triple CFa-styrene chain (between 7 and 



10 nm along the abscissa) images with decreasing height. 
Significantly, at Vs = —2.0 V, this region images 0.2 nm 
below the H:Si surface indicating depleted silicon state 
density beneath the molecules at the tip Fermi level. 

Results presented in Figures [l] and [2] and for other sin- 
gle chain CFs-styrene/OCHa-styrene heterostructures 
(not shown), show the presence of the chemical hetero- 
junction causes the interfacial OCHa-styrene to image 
with elevated height under low filled-state bias. The 
prominence (both in terms of height and lateral extent) 
of the filled-state interfacial feature compared with the 
absence of significant interfacial structure resolved in 
empty-state imaging in these structures, suggests an elec- 
tronic origin for the effect. As the low bias imaging con- 
ditions required to observe the filled-state interfacial fea- 
ture also result in decreased tip-sample separation, the 
possibility of tip-molecule interactions at the heterojunc- 
tion leading to altered molecular conformations which 
contribute to the observed interfacial structure cannot be 
ruled out. However, such tip-induced structural changes 
cannot account for the differing interfacial behaviour ob- 
served at low positive and negative substrate bias. The 
localized depletion in silicon filled-state density observed 
in Fig. [3] (typical of other side by side CFs-styrene 
structures studied) in response to the triple CFs-styrene 
chains also cannot be understood in terms of only confor- 
mational differences between otherwise non-interacting 
molecules. The theoretical modeling developed in the 
following Sections explores various factors which can ac- 
count for these observations. It will be seen that con- 
formational details in concert with electrostatics play a 
significant role in these phenomena. 



III. MODEL 

A. Formalism 

In order to carry out calculations of electronic trans- 
port through molecules bonded chemically to metal or 
semiconductor electrodes it is necessary to know the 
electronic structures of these systems. Ab initio den- 
sity functional calculations based on the Kohn-Sham lo- 
cal density approximation {LDAj^ and extensions such 
as generalized gradient approximations (GGA) are com- 
monly used for this purpose. However, the theoreti- 
cal foundations of this approach and the accuracy of 
the results that it yields for molecular transport calcu- 
lations are increasingly being questioned in the litera- 
ture at the present time:3 1° 11^12 13 14 15 16 n is 19 20 21 22] 

While such formalisms are appropriate (at least in prin- 
ciple) for calculating the total energy, the spatial dis- 
tribution of the electronic charge density and the elec- 
trostatic potential throughout i nhom ogeneous electronic 
systems in their ground states^^^^ the single-electron 
eigenenergies and wave functions that appear in them 
are somewhat artificial constructs tha t in mos t cases 
do not have rigorous physical meanings .I^Ii^IiSlll] Conse- 
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quently, they need not be good approximations to the 
energies and wave functions of the true electronic quasi- 
particles that determine the electronic transport prop- 
erties of molecular systems. Under what circumstances 
such ab initio calculations should produce acceptable 
results for electron transport despite this fundamental 
shortcoming, and how to obtain satisfactory results when 
they do not are important unresolved questions that 
are t he subject of much ongoing resear ch at the present 

^^Jj^(; l3ll0lliri2ll3ll4ll5ll6ll7ll8ll9l20t21l22l42l 

For molecules adsorbed on silicon the above deficien- 
cies of the LDA, GGA and their relatives manifest them- 
selves most obviously in that these approximations un- 
derestimate the band gap of silicon and yield incorrect 
values for the energy offset between the highest occupied 
molecular orbital (HOMO) of the molecule (or the rel- 
evant frontier orbital of the adsorbate) and the silicon 
valence band edge. The errors in these energy offsets ob- 
tained from the density functional calculations have re- 
cently been estimated for a few molecules to range from 
0.6 to 1.4eV.'^ Because of these and 

otheiHa! deficien- 
cies, the predictions of transport calculations based on 
ab initio density functional calculations of the electronic 
structure are unreliable for molecules on silicon; they are 
able to capture some observable phcnomcngp^ but are 
qualitatively incorrect for others.'^^'^ 

In the case of linear chains of styrend^ and methyl- 
styren molecules grown by self-assembly on a hydro- 
genated (001) silicon surface, the incorrect offsets be- 
tween the molecular HOMO levels and the silicon va- 
lence band edge given by the density functional calcu- 
lations result in such calculations yielding qualitatively 
incorrect STM images for these systems.^ In particu- 
lar, at low bias the ab initio calculations predict minima 
in the STM height profiles of the molecular chains over 
the centers of the molecules where maxima are observed 
experimentally.!^ These deficiencies of the ab initio cal- 
culations have been overcom^^H by developing a different 
electronic structure model based on extended Hiickel the- 
ory, a tight binding scheme from quantum chemistrjSSI^ 
that provides an approximate description of the elec- 
tronic structures of many molecules and has also been 
used successfully to explain the experimental current- 
voltage characteristics of a variety m olecular wires con- 
necting metal electroded^^'^^'^Slinilil and to model th e 
band structures of a variety of crystalline solid^^^'^^lSSI 

Extended Hiickel theory describes molecular systems 
in terms of a small set of Slater-type atomic orbitals 
{|0j)}, their overlaps Sij = {(f>i\(f>j) and a Hamiltonian 
matrix Hij = {(l)i\H\(j)j) . The diagonal Hamiltonian el- 
ements Hii — ei are taken to be the atomic orbital ion- 
ization energies and the non-diagonal elements Hij for 
i ^ j are expressed in terms of e;, e^, Sij and phe- 
nomenological parameters chosen for consistency with 
experimental molecular electronic structure data. As is 
described in detail in Ref. |38j the standard extended 
Hiickel theorj0Sl was modified so as to also provide an 
accurate description of the band structures of the silicon 



substrate and the tungsten STM tip. The energy offset 
between the molecular HOMO and the silicon valence 
band edge (an adjustable parameter in the theorjl^ 
was assigned a physically reasonable value for which our 
transport calculations reproduced correctly the character 
of the height profile along the styrene and methylstyrene 
molecular chains observed experimentally in STM im- 
ages at low bias, i.e., apparent height maxima over the 
centers of the molecules. The model proved to be re- 
markably successful, accounting not only for this low bias 
property of the STM images, but also for several counter- 
intuitive features of the experimental data, including the 
experimentally observed reversal in the contrast between 
the styrene and methyl-styrene molecular chains with in- 
creasing STM tip bias, the observed increase in the ap- 
parent height of the molecules at the ends of the molec- 
ular chains relative to those far from the ends with in- 
creasing bias and the observed disappearance of the cor- 
rugation of the STM height profile along the molecular 
chains with increasing bias .13^ 

The tight binding model based on extended Hiickel the- 
ory that is described in Ref. SSjis adopted in the present 
work but with an important modification: A limitation 
of extended Hiickel theory is that in it the atomic or- 
bital energy depends on the type of atom on which 
orbital i is located but is not influenced by the pres- 
ence of other atoms in the vicinity. This is a reasonable 
first approximation for the styrene and methylstyrene 
molecules considered in Ref. |3HI since those molecules 
do not contain strongly charged groups. However for the 
OCHa-styrene and CFa-styrene molecules that will be 
considered here, there is strong charge transfer between 
carbon atoms and fluorine and oxygen atoms that results 
in significant electrostatic fields that should modify the 
atomic orbital energies of surrounding atoms. These 
electrostatic fields are included in the present model as is 
explained below. The resulting variation in the electro- 
static potential from molecule to like molecule along the 
molecular chains has direct and striking effects on the ex- 
perimental STM height profiles presented in Section |TT] 
as will be discussed in detail in the Sections that follow. 

While, as has been discussed above and in Refs. 
I3ll0lllll2ll3ll4ll5ll6livil8llt)llj0l2ll22[ the use of den- 
sity functional theory at the level of LDA or GGA for 
calculating electronic quasiparticle properties has no fun- 
damental justification, this criticism does not apply to 
calculations of the ground state electronic charge density 
distribution and the electrostatic potentials that that have 
been shown by Hohenberg and KohrP^ to be functionals 
of the charge density. Therefore the use of ab initio den- 
sity functional calculations to estimate the electrostatic 
contributions to the atomic orbital energies that are 
due to charge transfer between different atoms (but are 
neglected in extended Hiickel theory), while still involv- 
ing approximations, is justified at the fundamental level. 
In the present work these electrostatic corrections to ex- 
tended Hiickel theory were included in our model in the 
following way: 
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An ab initio calculation was carried ouP^ of the elec- 
trostatic potential Wn at the nucleus of every atom n of 
an atomic cluster that included a molecular chain with 
a total of 20 OCHa-styrene and CFs-styrene molecules 
in the geometry that they take on the silicon surface, a 
few layers of nearby Si atoms, and the H atoms needed to 
passivate the dangling bonds on the surface of this silicon 
cluster. A similar ab initio calculation was carried out^* 
of the electrostatic potential C/„ at the nucleus of each 
of these atoms in the absence of all of the other atoms. 
Thus 

En = -e(W„ - Un) (1) 

is an estimate of the contribution to the electrostatic en- 
ergy of an electron in an atomic orbital on atom n that is 
due to the presence of all of the other other atoms in the 
system, calculated self-consistently from first principles. 
I.e., it is the electrostatic contribution to the atomic or- 
bital energy that is neglected in extended Hiickel theory, 
as discussed above. This contribution was included in the 
present model by making the substitution ei ei + En 
(for each atomic orbital i of every atom n) in the di- 
agonal elements Ha — of the extended Hiickcl-likc 
model Hamiltonian obtained for this system as described 
in Refi38> Because the orbital basis used is not orthog- 
onal (as in standard extended Hiickel theory), the non- 
diagonal matrix elements of the model Hamiltonian were 
also adjusted according to Hij — Hij + Sij{En + -Em)/2 
for orbitals i and?' on atoms n and m as required for 
gauge invariance.ESThe effect of the bias voltage applied 
experimentally between the STM tip and the silicon sub- 
strate on the Hamiltonian matrix H.^j was included in the 
present model as in RefI551 

Since in the above theoretical approach a highly 
compute-intensive ab initio calculation of the effects of 
charge transfer is carried out only once for any molec- 
ular chain being studied and a smaller basis set is used 
in the transport calculations that follow, in addition to 
correcting some fundamental deficiencies of density func- 
tional theory as discussed above, the present approach 
is able to treat much larger systems than is practical to 
study by ab initio methods alone. This was crucial for 
the present work where fairly long molecular chains on 
Si needed to be studied and the system is not periodic 
along the molecular chain. 

The electric current flowing between the STM tip and 
silicon substrate via the adsorbed molecules was evalu- 
ated as described in detail in Ref. [35] by solving the 
Lippmann-Schwinger equation, determining from the so- 
lution the Landauer transmission probability T(E, V) 
that depends on the electron energy E and appUed bias 
voltage V and evaluating the Landauer formula 

W = 7^y_ dETiE,V)if{E,fis)-fiE,fi^)) (2) 

where / is the current, f{E,fj,i) — 
1 / {exp[{E — fii) / kT] + 1) and /ii is the electrochem- 



ical potential of the source (i = s) or drain (z = d) 
electrode .Efil 

The chains of OCHa-styrene and CFa-styrene 
molecules and associated silicon clusters that needed to 
be studied theoretically to make comparison with the ex- 
periment included too many atoms for it to be practical 
to determine their relaxed structures using ab initio den- 
sity functional theory calculations. Therefore a molecular 
mechanics method was adopted: Chains of 40 molecules, 
half of the chain OCHa-styrene and the other half CF3- 
styrene, on a hydrogenated silicon cluster were relaxed 
using the Universal Force Field modeP^ starting from 
a variety of initial configurations, and many different 
metastable relaxed structures were found. It is reason- 
able to suppose that many of them as well as intermediate 
structures between them were being sampled thermally 
in the experiment which was carried out at room temper- 
ature. Thus the limited accuracy of the structures given 
by the molecular mechanics method was deemed to be 
adequate for the purpose of the present study. Since 
the phenomena of interest in the present work were ob- 
served experimentally near the junction of the OCH3- 
styrene and CFa-styrene chains, the relaxed molecular 
chain was truncated to a subset of 20 molecules surround- 
ing the junction and the calculations of the electrostatic 
potential and of electrical conduction between the silicon 
substrate and STM tip via the molecules were carried 
out for this truncated chain of molecules and underlying 
H-tcrminatcd Si cluster with a dimerized (001) surface 
without further relaxation. The STM tip in the present 
work was modelled as in Ref [38] as a 15 tungsten atom 
cluster with a (001) orientation and a single terminating 
atom, and coupled to an electron reservoir by many ideal 
leads. 



B. Energy Level Ordering 

In the present theoretical approach, the energy off- 
sets between the HOMO levels of the OCHa-styrene and 
CFa-styrene parts of the molecular chain are given by 
the extended Hiickel theory with the parameterization 
described in Ref. |35J modified so as to yield an accurate 
band structure for silicon and to include the ab initio 
electrostatic corrections En as discussed in Section |III A| 
The same applies to the molecular LUMO levels and to 
the variation of the atomic orbital energies from molecule 
to molecule along the molecular chain. 

The molecular HOMO and LUMO levels for both types 
of molecules considered here reside primarily on the ben- 
zene rings of the molecules. Thus, because of the elec- 
tron withdrawing (donating) nature of the CF3 (OCH3) 
group, the 0CH3-styrene HOMO is higher in energy 
than the CFg-styrene HOMO, ESI and the CFa-styrene 
LUMO is lower than the 0CH3-styrene LUMO. Thus 
the HOMO of the molecular chain as a whole is located 
on its 0CH3-styrene part (as has already been mentioned 
in Section |ll| and the LUMO of the molecular chain as 
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a whole is located on its CFa-styrene part. 

However, as in Ref. [38l the energy offset between the 
molecular HOMO of the molecular chain and the sili- 
con valence band edge is an adjustable parameter of the 
present theory whose value is not known accurately: As 
was discussed at the beginning of Section [HI A| this off- 
set is not given correctly by calculations based on stan- 
dard density functional theories; its accurate determina- 
tion (like the determination of the energy offsets between 
the molecular levels and substrate Fermi levels for other 
molecular nanowire^^^ is a difficult unsolved problem of 
molecular electronics. However, based on the discussion 
of our experimental results that follows it is plausible 
that the molecular HOMO is located below the silicon 
valence band edge and the molecular LUMO is located 
above the silicon conduction band edge: 

An important aspect of the experimental data pre- 
sented in Section |TT] is that even for the smallest values of 
\Vs\, the bias voltage voltage between the STM tip and 
substrate at which imaging of the molecular chains was 
feasible, \eVs\ was considerably larger than the band gap 
of silicon, for both positive and negative biasing of the 
silicon substrate relative to the STM tip. Given that the 
Fermi level of the STM tip at zero bias is located within 
the Si band gap, it therefore follows that the STM tip 
Fermi level was well below the top of the silicon valence 
band or well above the bottom of the silicon conduction 
band when the molecules were being imaged for negative 
and positive biasing of the substrate relative to the STM 
tip, respectively. A reasonable interpretation of this fact 
(that is consistent with all of the experimental data and 
with the results of our theoretical modeling) is that the 
HOMO of the molecular chain is located in energy be- 
low the silicon valence band edge, that the LUMO of 
the molecular chain is above the silicon conduction band 
edge, and that the experimental imaging of the molec- 
ular chains themselves (in the present study) was being 
carried out for STM tip Fermi energies near or above 
the LUMO or near or below the HOMO of the molecu- 
lar chain, conditions under which enhancement of STM 
currents due to resonant or near resonant transport via 
the molecules is to be expected. This will be the view 
adopted in the remainder of this article. 

Thus three regimes will be considered theoretically for 
both positive and negative substrate bias: 

1. Very low bias at which the STM tip Fermi level lies 
between the Si conduction band minimum (valence band 
maximum) and the molecular LUMO (HOMO) energy 
levels. 

2. Low bias at which the STM tip Fermi level is 
slightly above (below) the lowest (highest) molecular 
LUMO (HOMO) level. 

3. Higher bias at which the STM tip Fermi level is 
near the top (bottom) of the lowest (highest) molecular 
LUMO (HOMO) band. 

As will be seen below, each regime is predicted to have 
its own characteristic signature in STM imaging that can 
be compared with our experimental data. These signa- 



tures do not depend qualitatively on the precise value 
assumed for the offset between the silicon valence band 
edge and molecular HOMO level. Thus comparison be- 
tween theory and experiment allows us to draw conclu- 
sions regarding the specific regimes in which the data was 
being taken and the physical mechanisms underlying the 
phenomena that were observed. 



C. Prototypical Geometries of the Individual 
Molecules Bound to Si Dimers 



Different conformations have been proposed for styrene 
molecules on H terminated Si(lOO)^ '^^ and it is reason- 
able to suppose that different atomic geometries are pos- 
sible for the substituted styrenes on Si that are studied 
here. One possibility is that the C atom of the molecule 
that bonds to the C atom that bonds to the Si is located 
over the "trench" between two Si dimer rows as is shown 
in FigQa) for OCHs-styrene and in Figjij^b) for CF3- 
styrene. This conformation is similar to that assumed for 
chains of styrene molecules on Si in Refs 9 and 38 and will 
be referred to as "the T-tethered geometry" . In another 
conformation that has also been proposed for styrene on 
Si^^ the C atom that bonds to the C atom that bonds 
to the Si is located over the the Si dimer row to which 
the molecule bonds as in FigQc); this alternate confor- 
mation will be referred to as "the D-tethered geometry." 
In Figga) and (c) the C atoms of the OCH3 groups are 
located further over the "trench" than the O atoms are, 
i.e., the OCH3 group also has the T orientation. How- 
ever, the opposite (D) orientation of the OCH 3 group 
with the C oriented towards the Si dimer to which the 
molecule bonds as shown in FigQd) is also possible. 

The molecular conformations described above are ideal 
cases: Molecule-molecule interactions in the molecular 
chains and thermal motion at room temperature are ex- 
pected to result in many intermediate geometries with 
the molecular benzene ring and the CO bond of the 
OCH3 group not being coplanar with each other or with 
the silicon dimer to which the molecule bonds and the 
OCH3 and CF3 groups being rotated through different 
angles about the bonds that link them to their respec- 
tive molecules. As will be seen below, these deviations 
from the ideal geometries need to be taken into consid- 
eration when modeling the experimental STM images of 
the molecular heterowires. However it will often be con- 
venient in the discussions that follow to classify molecular 
structures according to the idealized structures such as 
those in Fig|4](T or D tethered molecule, T or D oriented 
OCH3 group) that they most closely resemble. 



9 



IV. THEORETICAL RESULTS FOR A 
REPRESENTATIVE MOLECULAR CHAIN 
GEOMETRY 

A. Structure 

In this Section we present theoretical results for elec- 
trical conduction between the STM tip and Si substrate 
through a particular 20-molecule CF3-styrene/OCH3- 
styrene chain on the dimerized H-terminated (001) silicon 
surface. The relaxed geometry of a part of this molecular 
chain surrounding the CFs-styrene/OCHa-styrene junc- 
tion is shown at the top of Fig. |5] Atoms are colored as 
in Fig|4] Each molecule bonds to a Si atom on the same 
edge of the same Si dimer row (the far edge of the dimer 
row closest to the viewer) through a single C-Si bond as 
in Fig|4] For the specific geometry shown in Fig. [5] the 
tethering of all of the molecules to the Si as well as the ori- 
entations of the OCH3 groups are of the T type defined 
in Section IIII CI Some further details of this structure 
should also be noted since they shall prove to be crucial 
in making the connection between the theoretical results 
that will follow and key aspects of our experimental data: 
Most of the molecules lean somewhat in the direction 
away from the heterojunction and also have swung out- 
wards somewhat (away from the heterojunction) about 
the axes of the tethering Si-C bonds, as one might ex- 
pect for molecules subject to a net weak mutual steric 
repulsion. Due to the repulsion between the F atoms on 
different CFa-styrene molecules, the CF3 groups are all 
rotated through different angles about the bonds between 
the CF3 groups and the molecules' benzene rings, and the 
benzene rings themselves are tilted somewhat differently 
from molecule to molecule. Because of this the detailed 
structure of CF3-styrene chain is more complex and less 
well ordered than that of the OCHa-styrene chain. 



B. Theoretical current profiles at positive and 
negative substrate bias 

Fig. |5] shows representative results for the calculated 
current / flowing between the tungsten STM tip and sil- 
icon substrate via the CF3-styrene/OCH3-styrene chain 
described in the preceding paragraph at some positive 
substrate biases (i.e., for empty molecular state imaging) 
vs. STM tip position along the chain at constant tip 
height. The corresponding results for negative substrate 
bias (filled molecular state imaging) are shown in Fig. [6] 
The black bullets (red diamonds) at the top of each plot 
indicate the lateral locations of the carbon atoms of the 
molecular CF3 groups (O atoms of the OCH3 groups). 

In both figures, curves L (black) are calculated current 
profiles along the molecular chain in the very low bias 



Regime 1 introduced in Section IIIB Here the STM tip 
Fermi level lies between the silicon conduction band edge 
and the lowest state of the molecular LUMO band in 
Fig. [5] and between the silicon valence band edge and 



the highest state of the molecular HOMO band in Fig. [6] 
The red curves E are in the low bias Regime 2 of Section 
|niB| where the STM tip Fermi level is slightly above 
(below) the lowest (highest) molecular LUMO (HOMO) 
level in Fig. [s] (Fig. |6]). The blue curves H are in the 
higher bias Regime 3 in which the STM tip Fermi level 
is near the top the lowest molecular LUMO band in Fig. 
|5] and just below the bottom of the highest molecular 
HOMO band in Fig. [6] 

The main features of the current plots in Figs |5] and [6] 
can be understood qualitatively by considering the pro- 
files along the molecular chain of the calculated electro- 
static shifts En of the atomic orbital energies defined by 
Eqn.([T|) that enter into the present model as discussed in 
Section [HI Al The relevant information is summarized in 
Fig. |7](a). There the red curve (Ph) shows the average 
of En over the six carbon atoms of the benzene ring of 
each molecule (the scale is on the left axis) vs. the lat- 
eral position of the center of the ring along the molecular 
chain. Since the HOMO and LUMO of both the CF3- 
styrene and OCHa-styrene are located principally on the 
benzene ring, the red curve is an approximate guide as 
to how the molecular HOMO and LUMO orbital ener- 
gies are affected by the electrostatic potentials due to 
charge transfer between the atoms of this system. The 
blue dotted vertical line separates the 0CH3-styrene part 
of the chain on the right from the CFa-styrene part on 
the left. The black curve labeled Si (for which the scale 
is on the right axis) shows _E„ for the Si atoms to which 
the molecules bond. 

These electrostatic shifts differ by several tenths of 
an eV between the 0CH3-styrene and the CF3"Styrene 
parts of the molecular chain due to the different signs 
and magnitudes of the charge transfer between the CF3 
and OCH3 substituent groups and the benzene rings of 
the respective molecules. They result in a further low- 
ering of the HOMO and LUMO of the CF3-styrene rel- 
ative to those of the 0CH3-styrene in addition to that 
predicted by extended Hiickel theory and also a very lo- 
cal band bending of the silicon conduction and valence 
bands in the immediate vicinity of the molecular chain; 
the Si band edges are lower in energy near the CF3- 
styrene chain than near the the 0CH3~styrene. The hor- 
izontal dotted line in Fig. |7](a) indicates En for surface 
Si atoms if all molecules are replaced by H atoms bound 
to the Si; the effect of such (passivating) H atoms on 
the orbital energies of the surface Si atoms is interme- 
diate between those of 0CH3-styrene and CF3-styrene 
molecules bound to the silicon. 

Because of the form of these electrostatic profiles, when 
the Si substrate is biased positively relative to the STM 
tip, the Fermi level of the metal tip crosses the Si con- 
duction band edge near the CFa-styrene chain and the 
LUMO levels of the CF3-styrene before doing so at the 
OCHs-styrene. Thus the CFa-styrene appears higher 
than the 0CH3-styrene for positive substrate bias. This 
is seen in the theoretical current plots in Fig. [5] and also 
experimentally in Fig. |2jc) although, as was discussed in 
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Section |lT] the contrast between the OCHs-styrene and 
CFs-styrene depends on the unknown microscopic de- 
tails of the STM tip. As will be discussed in Section [V| it 
also depends on the specific conformation of the molecu- 
lar chain. 

The red curve (Ph) in Fig. [7](a) has a sharp minimum 
at the third CFa-styrene molecule from the junction with 
the OCHa-styrene chain. Thus the molecular LUMO 
level is lowest at that molecule. Therefore as the (posi- 
tive) substrate bias increases, the Fermi level of the STM 
tip first crosses a molecular LUMO-derived state in the 
vicinity of that molecule. The prominent peak centered 
at the third CFa-styrene molecule in the red current pro- 
file E in Fig. |5]thus corresponds to the onset of electron 
injection from the STM tip into the molecular LUMO 
band. There is no such feature at this molecule in the 
black curve L in Fig. [s] (which is at lower bias) because 
in that case the STM tip Fermi level is still too far away 
from the molecular LUMO levels for resonant tunneling 
through those levels to be important; thus the larger elec- 
trostatic shift of the third CFa-styrene molecule's LUMO 
level does not play a large role here. When the bias is 
large enough for the STM tip's Fermi level to have risen 
to the top of the molecular CF3- styrene LUMO band all 
of the molecular LUMO orbitals throughout the CF3- 
styrene chain are transmitting resonantly so that the dif- 
ference between the third CFa-styrene molecule and the 
other CF3-styrene molecules has been largely erased in 
the blue curve H in Fig. [5] 

Notice that although the Si band edge is lowered lo- 
cally around the CFa-styrene molecules (see the black 
curve in Fig. |7] (a)) so that one may expect the Si to 
appear higher in that vicinity than elsewhere in empty 
state STM images, the black curve labelled Si in Fig. [7] 
(a) does not show a minimum at the third CFa-styrene 
molecule from the junction, confirming that it is the lo- 
cal electrostatic lowering of the molecular LUMO at that 
molecule and not the effect of the electric fields at the Si 
that is responsible for the prominent feature at the third 
CFa-styrene molecule in the red profile E in Fig. |5] 

At the bias voltages considered in Fig. [5] the tip Fermi 
level is well below the OCHa-styrene LUMO states so 
that resonant transport via the OCHa-styrene LUMO 
states is not important and thus there is little contrast be- 
tween different parts of the OCH3- styrene chain except 
possibly at the highest bias value shown for the OCH3- 
styrene molecule that is closest to the CFs-styrene chain. 

Notice also that in Fig. [5] the contrast between 
the CF3-styrene and 0CH3-styrene chains decreases 
markedly with decreasing magnitude of the applied bias. 
This is because the role of resonant transmission via the 
molecular levels in the CF3-styrene part of the chain is 
decreasing rapidly whereas no resonant transmission is 
involved in conduction via the 0CH3-styrene which is 
therefore less sensitive to the magnitude of the applied 
bias. This behavior is in qualitative agreement with the 
experimental data in Fig. [2]jc) and (d). 

The qualitative features of the theoretical current pro- 



files for negative substrate bias in Fig. [6] can also be un- 
derstood by considering the electrostatic energy profiles 
in Fig.|7](a): 

The current profile at very low bias (black curve L, 
Fig. |6| shows only weak contrast between different parts 
of the CF3-styrene chain and between different parts of 
the 0CH3-styrene chain because the tip Fermi level is 
between the silicon valence band edge and the molecular 
HOMO levels so that resonant tunneling via the molec- 
ular HOMO levels (that is strongly modulated by local 
extrema of the electrostatic potential along the chain of 
molecular benzene rings exhibited by the red curve in 
Fig. [7] (a)) is not occurring.'^ 

When the tip Fermi level falls below the highest molec- 
ular HOMO-derived level (which is located at and around 
the second 0CH3-styrene molecule from the junction 
where the electronic electrostatic potential energy [the 
red curve B in Fig. |7] (a)] has its maximum) the current 
develops a strong maximum there because of the onset 
of resonant conduction via the molecular HOMO in that 
vicinity. This is why the red curve E of Fig. [6] has a 
strong maximum at the second 0CH3-styrene molecule 
from the junction. 

With further increase of the bias voltage the tip Fermi 
level moves deeper into the molecular HOMO band and 
resonant transport via molecular states derived from the 
molecular HOMO becomes possible in other parts of the 
OCHs-styrene chain as well and thus the peak in the 
current profile of the OCHa-styrene chain becomes less 
pronounced, as in the blue curve H in Fig. [6) Notice 
also that with increasing bias the boundary between the 
OCHs-styrene and CF3-styrene chains becomes blurred 
with the higher current typical of the 0CH3-styrene 
chain extending to the nearest one or two CF3-styrene 
molecules, as in the feature labelled * in the blue curve 
H in Fig. |6l 

Another feature of Fig. |6]is the reversal in the contrast 
between the 0CH3-styrene and CF3~styrene chains be- 
tween low and intermediate bias: Curve L is higher for 
the CFa-styrene while curves E and H are higher for the 
OCHa-styrene. This is possible because the transport 
mechanisms for the OCHa-styrene in the low and inter- 
mediate bias regimes are different: The latter is domi- 
nated by resonant conduction via the molecular HOMO 
levels while the former is not. Whether such a contrast 
reversal with increasing bias occurs or not was found in 
the present theoretical study to depend on the structures 
of the molecular chains: For the geometry shown at the 
top of Fig. [5] ( see also Fig. |4]) the OCH3 groups lie espe- 
cially low relative to the CF3 groups which makes it pos- 
sible for the CF3-styrene molecules to appear higher than 
the 0CH3-styrene at low negative substrate bias. This 
occurs despite the fact that both the molecular HOMO 
levels and the local Si valence band edge are higher 
in energy on the 0CH3-styrene side of the heterowire. 
For some other geometries of the OCH3-styrene/CF3- 
styrene chain to be discussed in Section |V] the OCH3 
groups lie higher relative to the CF3 groups and no con- 
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trast reversal is found theoretically: The OCHa-styrene 
appears higher than the CFa-styrene even at low nega- 
tive substrate bias. In experiment, contrast reversal be- 
tween the 0CH3~styrene and CFs-styrene chains with 
changing bias was observed in some runs (as in Fig. [2| 
but not in others (as in Fig. [ij, most likely due to unde- 
termined differences between the STM tips involved. 

A significant difference between the calculated current 
profiles in Figs [S] and [6] and our experimental STM height 
profiles is that the latter do not show the quasiperiodic 
modulation with a period of roughly two molecular spac- 
ings seen in the calculated current profiles for the CF3- 
styrene part of the chain. This difference can be un- 
derstood as follows: This quasi-periodic modulation of 
the calculated CFa-styrene chain current profile is a ge- 
ometrical effect due to the differing orientations of the 
CF3 groups and also the associated differing tilts of the 
molecules along the CFa-styrene chain described near the 
end of Section |IV A| Therefore it is reasonable to expect 
this modulation not to be visible in STM images taken 
at room temperature where the CF3 groups are rotating 
rapidly about their axes and the local phase of the spa- 
tial current modulation is therefore fluctuating rapidly 
so that only an average over many structures with dif- 
ferent phases is observed on the time scale on which the 
STM scans are recorded, and thus a period equal to the 
molecular spacing is observed. 



C. The current peaks near the 
CFa—styrene/OCHa— styrene interface 

The prominent peaks in the red current profiles E in 
Figs |5] and [6] near the interface between the CFa-styrene 
and OCHa-styrene molecular chains have the following 
important characteristics: The peak for positive sub- 
strate bias in Fig. [s] is narrow, just one molecule wide 
and occurs on the CFa-styrene side of the junction. The 
peak at negative substrate bias in Fig. [6] is broader, ex- 
tending over a few molecules and occurs on the OCH3- 
styrene side of the junction. These properties of the in- 
terfacial peak at negative substrate bias are consistent 
with the experimental data presented in Section [llj Fur- 
thermore the interfacial peak at negative substrate bias 
was seen experimentally at low bias in all of the many 
CF3-styrene/OCH3-styrene heteroj unctions that we im- 
aged and persisted when the STM tip changed sponta- 
neously in the course of the imaging. By contrast, in 
empty state imaging (i.e., for positive substrate bias), 
no discernible height enhancement at the interface was 
usually observed. 

The above discussion in Section HV Bl identified the in- 
terfacial peaks in the calculated current profiles as be- 
ing due to a maximum or minimum of the electrostatic 
contribution to the molecular HOMO or LUMO energy 
along the molecular chain occurring at the location of the 
peak. However, in order to understand the physics un- 
derlying the differences between the predicted properties 



of the peaks at positive and negative substrate bias and 
the similarities and differences between the predictions 
and experiment it is necessary to examine how these fea- 
tures of the current profiles depend on the structure of 
the molecular chain. This will be explored in the next 
section. 



V. RELATIONSHIP BETWEEN STRUCTURE, 
ELECTROSTATICS AND CURRENT PROFILES 



Several significant aspects of the structures of the ex- 
perimentally realized molecular chains such as whether 
the molecules are tethered to the Si in the T or D orienta- 



tion (as defined in Section III C ) and whether the OCH3 
groups are T or D oriented, are not obvious from a direct 
inspection of the experimental STM images. It is also not 
clear theoretically which types of these structures should 
be favoured by the kinetics of the growth process or by 
energetics at room temperature. Thus it is desirable to 
study a variety of plausible relaxed molecular chain ge- 
ometries theoretically and to clarify how their structures 
relate to their electrostatic and current profiles. This is 
done in the present Section and the results presented here 
also elucidate the mechanisms underlying the phenomena 
described in Sections |TT] and |IV[ Potential energy profiles 
were calculated for many relaxed chain geometries and 
a few representative examples will be discussed below, 
along with possible experimental implications. 



A. Notation 

The important features of these molecular chain ge- 
ometries are summarized in Fig. |8] where projections of 
the positions of some key atoms of the molecular chains 
onto the Si (001) plane are shown: Each panel in Fig. [s] 
shows two rows of Si surface atoms belonging to different 
Si dimer rows. In each case the molecules are tethered to 
the lower row of Si atoms. The C atoms that bond to the 
Si atoms are orange. The C atoms that bond to the C 
atoms that bond to the Si atoms are violet. The C atoms 
of the CF3 and OCH3 groups are black. The O atoms are 
red and the F atoms are green. The C atoms belonging 
to the benzene rings that bond to O atoms or to C atoms 
not belonging to the benzene rings are white. The other 
other C atoms of the benzene rings, the H atoms and 
the other Si atoms are not shown for clarity. The dashed 
line shows the path of the tip atom of the STM for the 
theoretical STM current profiles for each structure that 
are presented in this paper. 

The structure considered in Section Hvl will henceforth 
be referred to as "Structure (a)" as per its label in Fig. [s] 
The other structures to be discussed here will be referred 
to similarly, according to their labelling in Fig. |8] Each 
of these structures consists of 20 molecules taken from 
the central region of a relaxed 40 molecule chain with- 
out further relaxation, as described at the end of Section 
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B. Structure (b) 

In Structure (b), as in Structure (a), the OCHa-styrene 
molecules are tethered to the Si in the T orientation. But 
the OCH3 groups are approximately reversed relative to 
Structure (a) and are now in the D orientation (defined in 
Section IIIC). In Structure (b) the tethering of the CF3- 
styrene molecules to the Si alternates between T and D 
and the orientation of the F atoms in the CF3 groups also 
alternates from molecule to molecule. However since this 
is a relaxed structure neither the OCHa-styrene chain nor 
the CFa-styrene chain has a truly periodic geometry. 



1. Electrostatic and current profiles 

The electrostatic potential energy profiles along the 
chain of benzene rings and along the row of Si atoms to 
which the molecules of Structure (b) bond are the red 
and black curves in Fig. |7jb), resp. Some representative 
current profiles for Structure (b) calculated at constant 
tip height along the dashed line in Fig. ^h) are shown in 
Fig. |9] together with a side view of a part of this molec- 
ular chain. The qualitative behavior of the calculated 
current profiles in Fig. [9] can be understood by consider- 
ing the features of the potential energy profile along the 
chain of benzene rings in Fig. [7|b) an d applying reason- 
ing analogous to that in Section ]lVB[ The current profile 
under a moderately low positive substrate bias for which 
the STM tip Fermi level is just above the lowest level of 
the molecular LUMO band (the violet curve E') has a 
strong maximum one molecule wide at the CFa-styrene 
molecule next to the OCHa-styrene chain where the ben- 
zene ring potential energy profile has its sharp minimum. 
The current profile under a moderately low negative sub- 
strate bias for which the STM tip Fermi level is just below 
the highest level of the molecular HOMO band (the red 
curve E) rises gradually to a broad maximum at the end 
of the OCHs-styrene chain remote from the CFa-styrene 
where the benzene ring potential energy profile has its 
broad maximum. Unlike in Figj6] in this case there is no 
contrast reversal between the OCHa-styrene and CF3- 
styrene chains with increasing negative substrate bias; 
the OCHa^styrene is higher in both the very low bias 
profile L and the moderately low bias profile E. This is 
because the OCH3 group is located higher relative to the 
CF3 group in Structure (b) than in Structure (a). 



2. Significance of the results 

The two main differences between the electrostatic pro- 
files for Structure (b) and those for Structure (a) (shown 
in Fig. |7|^b) and Fig. |7|^a)) are: 



1. The minimum of the benzene ring electronic po- 
tential energy profile (while remaining very sharp) has 
shifted from the third CF3-styrene molecule from the 
heteroj unction in Fig. [7](a) to the CF3~styr ene molecule 
next to the heteroj unction in Fig. [7]^b). 

2. The maximum of the benzene ring potential energy 
profile while remaining broad has shifted from the vicin- 
ity of the heterojunction (on the 0CH3-styrene side) in 
Fig. [7]^ a) to the vicinity of the end of the 0CH3-styrene 
chain that is remote from the heterojunction in in Fig. 

Systematic studies of various structures and their po- 
tential energy profiles showed the reversal of the slope 
of the potential energy profile along most the chain of 
benzene rings of the 0CH3-styrene chain from Structure 
(a) to Structure (b) in Fig. |7] should be attributed not 
to the OCH3 groups being switched from the T to the 
D orientation per se, but primarily to the fact that the 
projection onto the axis of the molecular chain of the 
dipole moment associated internal charge structures of 
the OCH3 groups (negative O and positive CH3) has re- 
versed direction from Structure (a) to Structure (b) as 
is evident from the geometries of the two structures in 
Fig. [8] This reversal of the charge structure along the 
molecular axis should lead to a reversal of the slope of 
the potential energy profile and this is what is in fact 
seen in Fig. |7] 

In this way the change in the geometrical structure of 
the OCHa-styrene chain has resulted in the shift of the 
potential energy maximum (and STM tip current max- 
imum at moderately low negative substrate bias) from 
the vicinity of the heterojunction to the vicinity of the 
far end of the OCHs-styrene chain. 

An important point to note is that because of the ge- 
ometrical structure of the OCHa-styrene molecules and 
the way in which they bond to the silicon the internal 
dipole moment of the OCH3 group has a strong compo- 
nent parallel to the silicon surface and the OCH3 groups 
of all of the 0CH3-styrene molecules can be aligned in 
such a way that the projections of the OCH3 groups of 
all of these molecules onto the axis of the molecular chain 
have the same sign. This results in a gradual buildup of 
the electrostatic potential along the 0CH3-styrene chain 
and a broad electrostatic potential energy peak near one 
of its ends as in Fig. [7](a) or (b), and consequently a 
broad current peak under moderately low negative sub- 
strate bias, consistent with the experimental results pre- 
sented in Section |IT1 

By contrast, because the geometrical orientation of the 
CF3 group (see Fig. |4]) and the arrangement of charges 
within it are different than for the OCHa group, and be- 
cause of the repulsion between the F atoms on different 
molecules (see the end of Section IV A I , such a well or- 
dered arrangement of dipole projections onto the axis of 
the molecular chain does not occur for the CFa-styrene 
chain. Thus the electrostatic potential energy minimum 
on the CF3-styrene side of the heterojunction is a much 
more local phenomenon. Its location is determined by 
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the orientations of a few molecules surrounding it and by 
the fringing electrostatic field of the OCHa-styrene chain 
that extends across the heterojunction into the CF3- 
styrene chain. Thus the electrostatic potential energy 
minimum on the CFs-styrene side of the heterojunction 
is very narrow and its location is sensitive to the details 
of the local molecular geometry. It is therefore reason- 
able to expect any profile peak associated with it in STM 
imaging at room temperature (at positive substrate bias) , 
if present at all, to be more strongly impacted by ther- 
mal fluctuations in the molecular geometry and there- 
fore less readily observed than the corresponding peak at 
negative substrate bias on the OCHs-styrene side of the 
molecular chain. This is consistent with the differences 
between the experimental results obtained under positive 
and negative substrate bias that are presented in Section 

im 

Experimentally, height enhancement was observed in 
the STM images towards the end of the OCHa-styrene 
chain that is far from the junction with the CFa-styrene 
at low negative substrate bias as can be seen, for ex- 
ample, in Fig. [ijd). Because in this work the growth of 
the 0CH3~styrene chain followed the growth of the CF3~ 
styrene, a charged dangling bond is normally expected to 
be present at the end of the OCHa-styrene chain (remote 
from the OCHa-styrene/CFa-styrene junction) where 
the growth of the molecular chain terminated, as in pre- 
vious experimental work on styrene chains on silicon.!^ 
The electrostatic potential due to this charged dangling 
bond gives rise to enhancement of the height profile of 
the molecular chain in its vicinity in STM experiments 
as is discussed in Refl6^ If alignment of OCHa-styrene 
molecular dipoles as in Structure (b) was occurring in the 
present experiments, the associated electrostatic profiles 
such as in Fig. ^h) would also contribute to enhanced 
current near the end of the molecular chain as in curve 
E in Fig. |9] and hence to the enhanced height profile 
there in constant current STM experiments as in Fig. 
[ijd) . Determining experimentally whether such a molec- 
ular dipolc contribution to the height profile in vicinity of 
the dangling bond is present or not is in principle possible 
by passivating the dangling bond through the addition 
of a H atom and observing whether any height enhance- 
ment remains afterwards near the end of OCHa-styrene 
chain. Carrying out such a test was beyond the scope 
of the present experimental work. However, the fact 
that height enhancement was consistently observed ex- 
perimentally near the OCHa-styrene/CFa-styrene het- 
erojunction at low negative substrate bias suggests that 
OCHa-styrene chains with dipole fields similar to that 
arising from Structure (a) rather than Structure (b) in 
Fig. 8 played the dominant role in these experiments. 



C. Structure (c) 

Structure (c) consists of 5 CFa-styrene and 15 OCHa- 
styrene molecules. The geometry of the CFa-styrene 



molecules and of the OCHa-styrene molecule that is next 
to the CFa-styrene is similar to that of the corresponding 
molecules in Structure (b). However, as can be seen in 
Fig. [sjc), the other OCHa-styrene molecules in Structure 
(c) have a different geometry: they are tethered to the 
Si in the D orientation and their OCHa groups are in the 
T orientation. Between the 8th and 9th OCHa-styrene 
the molecules from the junction with the CFa-styrene 
there is another abrupt change in the molecular geome- 
try: Although the molecules to the right and left of this 
dislocation are both D-tethered and the OCHa's are T 
oriented, the molecules to the right of the fault tilt more 
to the right and their OCHa groups are oriented nearly 
perpendicularly to the axis of the molecular chain. 

This is reflected in the electrostatic potential energy 
profile along the chain that is plotted in Fig. [tJ^c): Since 
the molecular structure at the junction is similar to that 
for Structure (b) the potential profile on the benzene 
rings is also similar there, i.e., there is a sharp poten- 
tial energy minimum at the CFa-styrene molecule at the 
junction. Since for OCHa-styrene molecules 2-8 from 
the junction the orientation of the O and two C atoms 
bonded to it is similar to that for the OCHa-styrene 
molecules in Structure (a) the orientation of the rele- 
vant dipoles on that part of the chain is also similar and 
thus the potential energy profile is also qualitatively sim- 
ilar: There is a broad potential energy maximum on the 
OCHa-styrene benzene rings peaked near the junction 
with the CFa-styrene. However, the different orienta- 
tion of the charged groups of the of the 9th- 15th OCHa- 
styrene molecules from the junction results in a reversal 
of the slope of the potential profile and the electronic po- 
tential energy on the OCHa-styrene benzene rings rises to 
another maximum near the right hand end of the molec- 
ular chain. 

The calculated STM current profiles for Structure (c) 
shown in Fig. [TO]are again consistent with the potential 
energy profile along the chain of benzene rings: At mod- 
erately low positive substrate bias (violet curve E') the 
current profile has a sharp maximum one molecule wide 
at the CFa-styrene molecule where the potential energy 
minimum along the molecular chain is located. At mod- 
erately low negative substrate bias (red curve E) there is 
a broad current maximum peaked on the OCHa-styrene 
side of the junction where the potential energy has one 
of its maxima. From there the current falls off to a min- 
imum near the middle of the OCHa-styrene chain where 
the potential energy curve also has its minimum. To 
the right of this minimum both the potential energy and 
current rise again but the current does not rise as high 
as near the junction and then begins to fall off again. 
This difference is due to the fact that for the red curve E 
the electron transmission through the molecules has not 
quite reached resonance for the OCHa-styrene HOMOs 
at the right end of the OCHa-styrene chain (although 
it has done this for the OCHa-styrene near the junc- 
tion) and also because the OCHa-styrene molecules are 
tilted more at the right end of the chain than at the left. 
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This geometrical difference also shows up as a step in the 



other current profiles in Fig. 10 at the dislocation in the 
OCHa-styrene chain. 

For Structure (c) there is again no reversal of the 
contrast between OCHs-styrene and CFs-styrene chains 
with increasing bias. 



D. Structure (d) 

In Structure (d) the geometry of the CFs-styrene chain 
is broadly similar to those of Structures (b) and (c). 
However, in the OCHa-styrene chain the tethering to 
the Si alternates between D and T. The orientation of 
the OCH3 groups is all T and most of the relevant OC 
dipoles are oriented roughly as for Structure (a). Thus 
one should expect a potential energy profile along the 
chain of 0CH3-styrene benzene rings somewhat simi- 
lar to that for Structure (a). Indeed the red curve in 
Fig. |7f^d) does have its maximum in the OCHa-styrene 
chain near the junction with the CFa-styrene and it de- 
clines overall to the right, but the decline is modulated 
by strong oscillations due to the different tethering of 
alternate OCHs-styrene molecules to the Si. The poten- 
tial energy profile along the chain of benzene rings on the 
CFs-styrene chain of Structure (d) is quite different than 
that for Structures (b) and (c) : There is no potential en- 
ergy minimum on the CFa-styrene next to the junction 
but instead a shallow minimum near the middle of the 
CFs-styrene chain. 

The current profiles for Structure (d) are shown in Fig. 
|ll| for moderately low positive substrate bias (violet curve 
E') and moderately low negative substrate bias (red curve 
E). Again the maxima of these profiles are near the lowest 
and highest points of the potential energy profile along 
the chain of benzene rings, respectively. But in this case 
the potential energy peak in the OCHa-styrene chain is 
less prominent and narrower than in the the previously 
discussed cases due to the less uniform structure of the 
array of dipoles. The broad potential energy minimum 
in the CFa-styrene results in a broader current peak at 
positive substrate bias than in the cases discussed earlier. 

This illustrates further the strong sensitivity of the 
CFs-styrene chain's potential minima to the precise de- 
tails of the conformation of a single molecule or a few 
molecules: In Structure (d) the OCHa-styrene molecule 
at the junction of the chains is D-tethered and the OCH3 
group is T oriented while the reverse is true for Structures 
(b) and (c) while the CF3 groups are rotated slightly dif- 
ferently and CFa-styrene at the junction tilts somewhat 
differently. 



E. Structure (e) 

Structure (e) is similar to Structure (a) except that 
now the orientations of the OCII3 groups alternate be- 
tween T and D and the CF3 groups have rotated about 



their axes through substantial angles. Because of the 
alternating orientations of the OCH3 groups the elec- 
tric fields due to the different 0CH3^tyrene molecules 
are not aligned and produce a potential energy profile in 
Fig. [Tje) with a very broad maximum near the end of 
the chain remote from the junction. The different orien- 
tations of the CF3 groups than in Structure (a) result in 
a sharp potential energy minimum at the second CF3- 
styrene molecule from the junction. 

The current profiles for Structure (e) are shown in Fig. 



12 for moderately low positive substrate bias (violet curve 
W) and moderately low negative substrate bias (red curve 
E). As expected for such a potential energy profile we 
see a strong narrow current peak on the second CF3- 
styrene molecule from the junction. The broad potential 
energy maximum on the OCHs-styrene chain produces a 
broad current maximum peaked towards the end of the 
0CIl3-styrene chain that is remote from the junction at 
intermediate negative substrate bias. 

For both Structure (e) and Structure (d) the alternat- 
ing geometries along the OCHs-styrene chains result in 
a strong modulation of the current with a period of 2 
molecular spacings. Since this modulation is not seen 
experimentally, either these structures are not realized 
in the experiments or the experiments are averaging over 
geometries for which the modulation occurs with differ- 
ent phases due to thermal fluctuations of the molecular 
geometries at room temperature. 



VI. ELECTROSTATIC VS. ELECTRONIC 
MOLECULE-MOLECULE COUPLING 

For comparison, similar calculations to those described 
above were also carried out for Structure (a) 

(i) omitting the environmental electrostatic shifts of 
atomic orbital energies En defined by Eqn. [T] No en- 
hancement of the current near the heterojunction (such 
as m curves E in Fig. |5]and[6| was found in those calcu- 
lations. 

(ii) Transport calculations were also carried out includ- 
ing the electrostatic shifts £'„ but with all Hamiltonian 
matrix elements and basis function overlaps responsible 
for direct electronic hopping from molecule to molecule 
switched off. These calculations produced very similar 
current profiles along the molecular chain to those in Fig. 
[5] and [6] including the interface current peaks and even 
the feature labelled ★ in the blue curve H in Fig. |6] How- 
ever these features appeared at somewhat different values 
of the bias voltages because the molecular HOMO and 
LUMO bands are narrower (as expected) when the inter- 
molecule Hamiltonian matrix elements and basis state 
overlaps are turned off. 

These results confirm that the enhanced current fea- 
tures near the interface in Fig. [s] and [6] result from 
electrostatic fields established by the polar molecules. 
These results therefore suggest an electrostatic origin for 
the fiUed-state interfacial enhancement resolved exper- 
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imentally in Fig. [T] and [2] at low bias. (Significantly, 
no such interfacial current enhancement was observed 
experimentalljEH in styrene/methylstyrene heterostruc- 
tures on Si where molecular electric fields are expected 
to be much weaker because of the weaker charge transfer 
within those molecules.) 

Furthermore, although electronic hopping from 
molecule to molecule may be occurring quite efficiently 
in substituted styrene chains on Si (100) as is discussed 
in Ref. |38l these calculations demonstrate the hopping 
not to be necessary for the occurrence of the interfacial 
features reported in the present work. 



VII. MODELING SINGLE-TRIPLE ROW 
CF3-STYRENE HETEROSTRUCTURES AND 
THEIR INFLUENCE ON CONDUCTION 
THROUGH THE UNDERLYING SILICON 

The main features of the calculated current profiles 
of the molecular heterostructures modeled theoretically 
in the preceding sections derive directly from the elec- 
tronic structures of the molecules themselves modulated 
by the electrostatic fields due to other nearby molecules. 
This however is not the case for the most striking fea- 
ture of the experimental data shown Fig. [Sj^d) and (e) 
where at low negative substrate bias the electric current 
passing through a triple row of CFa-styrene molecules 
is strongly depressed relative to that passing through the 
nearby single row of CFa-styrene molecules. Our theoret- 
ical findings presented below suggest that the influence of 
the adsorbed CFa-styrene molecules on the electrostatic 
potential in the underlying silicon may be responsible for 
this phenomenon. 

Representative results of our calculations are presented 
in Fig. 13 The molecular geometry studied is shown 
schematically in Fig 13 ^a). In the model considered the 
long CFs-styrene line is located between the two short 
ones to minimise sensitivity to the Si cluster edges. All 



of the molecules are assumed to be tethered to the silicon 
in the T orientation defined in Section [ill Gl and shown in 
Fig. Qb). Fig|T3|^b) shows the simulated constant height 
current along the long (central) CFa-styrene line. At low 
negative substrate bias (for which the STM tip Fermi 
level is close to the silicon valence band edge) , the calcu- 
lated current (curve V) drops by a factor of ~26 from the 
tallest current peak near the left hand end of the single 
CFa-styrene line to the tallest peak near the right hand 
end of the central line of the triple CFs-styrene structure. 
The origin of this effect is seen in the solid black curve 
labeled "Si" in Fig 13 'c): Dipole fields associated with 



the CFs-styrene molecules lower the Si orbital energies 
below the triple CFa-styrene by ~0.2eV more than under 
the single file CFa-styrene. At low bias this reduces the 
silicon electronic density of states near the silicon surface 
at the Fermi level of the STM tip more under the triple 
CF3-styrene line than under the single CFs-styrene line, 
resulting in the weaker tip current through the former. 



As the magnitude of the bias increases eventually the 
tip Fermi level is lowered further by an amount substan- 
tially exceeding the differential electrostatic energy shift 
between the Si orbitals under the triple and single CF3- 
styrene lines and thus the contrast between the STM tip 
currents through the triple and single CFa-styrene lines 
becomes weaker as is seen in curve L of Fig(l3];b) for 
which the tip Fermi level is 0.5eV lower than for curve 
V. 

This behavior found in the simulations (lower current 
through the triple CFs-styrene line than through the sin- 
gle CFs-styrene line at low bias and the reduction in this 
contrast with increasing magnitude of the bias) is quali- 
tatively similar to that seen experimentally in Fig. ^e). 
However, quantitatively the suppression of the current 
through the triple CFa-styrene line at low bias is more 
pronounced in the experiment (Fig. [3|e)) than in the 
theoretical curve V of Fig|T3|b): A 0.2-0.3nm difference 
in tip height at constant current (as between the triple 
and single CFs-styrene lines in Fig. |3je) at low bias) 
normally corresponds to approximately a 2-3 orders of 
magnitude change in current at constant tip height in 
STM experiments. 

Part of this difference between the simulation and ex- 
periment may be due to the unknown microscopic details 
of the STM tip, which are known to influence the height 
contrast between CFs-styrene and OCHs-styrene lines as 
is discussed in Section [iT] and can be seen in Fig(2] How- 
ever, the limited size of the model system for which the 
present simulations could be carried ouif^^l is clearly re- 
sponsible for at least a part of the difference between the 
experimental and theoretical findings: While the down- 
ward electrostatic shift of the Si orbital energies (as in- 
dicated by the solid black curve in Fig 13 'c)) is clearly 
larger under the triple than the single CFs-styrene line, 
it is evident that the 8 molecule single CFa-styrene line 
modeled is not long enough for the electrostatic poten- 
tial in the silicon under the single CFa-styrene line to 
reach a plateau with increasing distance from the junc- 
tion with the triple CFs-styrene line, i.e., the electro- 
static influence of the triple row of molecular dipoles of 
the triple CFa-styrene line is clearly being strongly felt 
under much of the single CF3-styrene line as well, re- 
sulting in diminished contrast in the calculated current 
between the single and triple CFs-styrene lines relative to 
that which may be expected in larger model systems. In 
this regard it is instructive to compare the electrostatic 
profile in the silicon under the single-triple CFs-styrene 
heterostructure with that for a somewhat longer single 
file CF3-styrene/OCH3-styrene heterostructure shown by 
the dashed black curve in Fig 13 'c): For the latter het- 
erostructure the potential profile in the silicon under the 
CF3-styrene line (which is on the left) does show a pro- 
nounced plateau where the value of the electrostatic po- 
tential is close to that for the single-triple CFs-styrene 
heterostructure under the end of the single CFs-styrene 
line that is furthest from the triple CF3-styrene line. 

In comparing theory with experiment it is noted that 
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there is also contrast between the tip currents calculated 
for tip positions over the three molecular rows within 
the triple row structure: The calculated current for the 
tip over the leftmost molecular row in Fig 13 'a) for the 
same tip height and bias as for plot V of Fig 13 'b) is 
shown by the dashed curve at the bottom of Fig 13 'b): 
The current over the left molecular row is smaller than 
over the center row of the triple structure (plot V) . The 
calculated current over the rightmost row of Fig|T3|(a) 
(not shown) is intermediate between that for the left and 
center rows. These results suggest that stronger con- 
trast between the single CFs-styrene line and parts of 
the triple CFa-styrene than is seen in plot V of Figp^b) 
may be possible theoretically. However, calculations for 
larger silicon clusters than were feasible in this study are 
required to eliminate possible cluster edge effects and de- 
velop a better understanding of this issue. 



If the magnitude of the negative substrate bias is in- 
creased still further, well beyond that for current plot L 
in Figp^b), the STM current through the single CF3- 
styrene line is predicted by the present simulation to be- 
gin to increase more rapidly than that through the central 
line of the triple CFs-styrene structure. This is seen in 
curve H of Fig |13[ b) where the current through the single 
CFs-styrene line is again much larger than that through 
the triple CFa-styrene line. This effect can be under- 



stood by considering the red curve (Ph) in Fig 13 c): 
This shows the average of En over the six carbon atoms of 
the benzene ring of each molecule of the long CFs-styrene 
row. This is higher for the single CFa-styrene line than 
for the triple CFa-styrene line. Therefore the molecu- 
lar HOMO levels for the molecules in the single CF3- 
styrene line are higher in energy than those for molecules 
in the center row of the triple CF3-styrene line. Con- 
sequently the molecular HOMO levels for the molecules 
in the single CFs-styrene line are crossed by the Fermi 
level of the STM tip at lower magnitudes of the nega- 
tive substrate bias than for molecules in the triple CF3- 
styrene line which results in the renewed enhancement 
of the current through the single CFa-styrene line rel- 
ative to the triple CFa-styrene line seen in curve H of 
Figp^b). No recovery of stronger contrast between the 
single and triple CFa-styrene lines at stronger negative 
substrate bias is seen experimentally in Fig. indicat- 
ing that even for the strongest negative substrate bias val- 
ues realized in these experiments the STM tip Fermi level 
remains well above the CFs-styrene HOMO levels even 
for the molecules in the single CFs-styrene line. This is 
consistent with our experimental findings for the CF3- 
styrene/OCHa-styrene heterostructures where no accel- 
eration of the increase of the apparent height of the CF3- 
styrene with bias at strong negative substrate bias (that 
would be indicative of the STM tip Fermi level crossing 
the CFa-styrene HOMO) was observed. 



VIII. DISCUSSION 

The present experimental and theoretical work has 
begun to explore how electric fields emanating from 
molecules influence electrical conduction through other 
molecules in their close vicinity and in the underlying 
substrate. 

In this work we have employed density functional 
calculations only to estimate the electrostatic potential 
throughout the system in its electronic ground state, an 
application of density functional theory that is believed 
to be soundly based at the fundamental level. The phe- 
nomena that we report on here can be understood qual- 
itatively by considering just the results of these electro- 
static calculations as summarized in Fig. [T] and 13 c). 
Transport calculations are however necessary to trans- 
late the electrostatics into current profiles that can be 
compared with the experimental STM data. As is dis- 
cussed in Section [III A| transport calculations that incor- 
porate electronic structures obtained entirely from den- 
sity functional theory, although popular, are not soundly 
based and are often misleading, especially for molecules 
on semiconductors. We therefore base our transport cal- 
culations instead on electronic structures obtained from 
semi-empirical tight binding models modified to include 
the density functional theory-based ground state electro- 
static potentials. 

Our theoretical modeling indicates the observed low 
bias filled-state current enhancement in the interfacial 
0CH3-styrene molecules to be due to the collective effect 
of the electric dipolc fields generated by the OCH3 groups 
of the OCHs-styrene chain when the OCH3 groups are 
aligned preferentially so that their carbon atoms are fur- 
ther from the heteroj unction than their oxygen atoms, 
as they are, for example, in Fig. |8]^a). It is plausible 
that structures of this sort are prevalent in these sys- 
tems since some of them should be favored for molecules 
subject to a weak net mutual steric repulsion as is dis- 
cussed Section [IV A[ However the resolution of our exper- 
imental STM data which was taken at room temperature 
is not sufficient to determine independently whether or 
not the molecular chains have structural order of this 
sort. Furthermore it is apparent from a comparison of 
our experimental data and theoretical results that the 
CF3 groups in the CFa-styrene molecules must be spin- 
ning rapidly enough on their axes that only an average 
over many conformations of these molecules is being ob- 
served in the STM images. Similarly it seems possible 
that many conformations of the 0CH3-styrene molecules 
are being averaged over rapidly in the STM imaging and 
that the above alignment of the OCH3 groups is simply 
more commonly present than other structures, such as 
that in Fig. [8|b) for which the apparent height enhance- 
ment should be near the end of the OCHs-styrene chain 
that is remote from the junction. 

The behavior of the imaging height enhancement that 
we observe experimentally with increasing negative sub- 
strate bias is also consistent with the predictions of our 
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theoretical model of this phenomenon. 

Under positive substrate bias simulations indicate that 
height enhancement near the junction should also oc- 
cur but should in this case be confined to only a single 
molecule. The location of the predicted height enhance- 
ment under positive substrate bias is sensitive to the de- 
tails of the local molecular geometry in its vicinity that 
are expected to fluctuate rapidly at room temperature. 
Thus definitive observation of this effect will most likely 
require cryogenic experimental work. 

The work reported here has succeeded in identify- 
ing transport phenomena occurring at molecular length 
scales that can reasonably be attributed to electric fields 
emanating from polar molecules impinging on other 
nearby molecules and on the underlying substrate, and in 
developing some insights into the detailed physical mech- 
anisms that may be involved. As has been discussed in 
Section|l]these phenomena may ultimately find device ap- 
plications, such as molecular switches in which the cur- 
rent through a molecule is switched by conformational 
changes in another nearby molecule and/or devices in 
which conduction through a semiconductor substrate is 



controlled with the help of molecules at its surface. 

However further experimental and theoretical work is 
needed in order to better understand and ultimately con- 
trol such phenomena. 
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FIG. 1: Color online. Constant current filled-state STM images showing the growth of two CFa-styrene/OCHa-styrene 
heterowires. (a) Vs = —3.0V. H:Si(100) surface following a lOL exposure of CFa-styrene. Arrows indicate positions of 
chemically reactive dangling bonds at the ends of the CFs-styrene line segments, (b) Following a lOL exposure of OCH3- 
styrene, the two CFa-styrene line segments in (a) have been extended to form two CFa-styrene/OCHa-styrene heterowires 
('1' and '2'). At K = -3.0V, the OCHa images higher (brighter) than the CFa-styrene segments, (c) 14 = -1.8V. The 
OCHa-styrene near the molecular heterojunctions in heterowires 1 and 2 image with enhanced height. Molecules at the end 
of the OCHa-styrene line segment in heterowire 1 also image with enhanced height in response to the terminal dangling bond. 
Inset: Close-up of heterowire 1. Molecules are bound to the right-hand side of the dimer row indicated by the red dots, (d) 
Constant-current topographic cross-sections extracted from bias dependent imaging of heterowire 1 along the trench to the 
right of the attachment dimers in (c). The height envelope for the heterostructure extends between 1 nm and 9.5 nm along the 
abscissa. The height maxima associated with the terminal dangling bond, and the molecular heterojunction are at 2.3 nm, and 
6.4 nm, respectively. At elevated bias (-3 V), the OCHa-styrene images with approximately uniform height from beyond the 
terminal dangling bond to the heterojunction. As the bias decreases in magnitude, the OCHa-styrene images with decreased 
height as the molecular tt states drop below the tip Fermi-level. At -1.8 V, localised height enhancement in the OCHa-styrene 
due to the terminal dangling bond and also near the molecular heterojunction (black arrow) is most evident. Image areas 
(a)-(c): 26nm x 26nm. Inset area: 8.5nm x8.5 nm. Tunnel current: 40 pA. 
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FIG. 2: Constant current fiUed-state (a)-(b) and empty-state (c)-(f) STM imaging of a 3 heterowire cluster, (a) Vs = —2.6V. 
OCHs-styrene line segments image higher (brighter) than CFs-styrene line segments, (b) Vs = —2.0V. At low magnitude 
bias OCHs-styrene can image lower than CFa-styrene (tip dependent). Enhanced molecular conductance throughout the 
interfacial OCHs-styrene remains evident, (c) Vs = -I-2.6V. Consistent with the greater electron affinity of CFs-styrene, 
CFa-styrene images with increased height relative to OCHs-styrene. (d) Vs = -I-2.0V. At reduced bias, height contrast 
between the CFs-styrene and OCHs-styrene line segments decreases, and molecules image with similar corrugation, (e) and 
(f) Vs = +2.0V. Depending on tip structure, OCHs-styrene can image above CFs-styrene and with varying (tip dependent) 
corrugation. Greyscale: nm (black), and 0.59nm, 0.44nm, 0.28nm, O.lSnm, 0.39nm, O.lSnm (white), (a) to (f), respectively. 
Image areas (a)-(f): 26nmx26nm. Tunnel current: 40pA 
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FIG. 3: Color online. Constant current filled-state STM images showing the growth of a (single-triple CF3-styrene)/OCH3- 
styrene heterostructure. (a) Vs ~ —3.0V. H:Si(100) surface following a 10 L exposure of CFa-styrene. White arrow indicates 
the position of the chemically reactive dangling bond at the end of the long CFa-styrene segment. The white asterisk shows 
a short double CFs-styrene line segment beside the longer single CFa-styrene chain, (b) Vs = +2.0V. Following a 10 L 
exposure of OCHs-styrene, the long CFs-stryene chain has been extended by ~ 7 OCHs-styrene molecules, (c) Va = —3.0V. 
OCHs-styrene images above CFa-styrene. Single and triple CFs-styrene segments image with similar height, (d) Vs = —2.0V. 
Single chains of OCHs-styrene and CFs-styrene continue to image above the H:silicon surface (brighter). Triple chains of 
CFa-styrene image below the H:silicon surface (black), (e) Constant current topographic cross-sections (0.4nm wide) extracted 
from bias dependent imaging of the CFs-styrene/OCHs-styrene heterowire along the trench to the right of the attachment 
dimers. Heights are given relative to the H:silicon surface (height = nm). Image areas (a)-(d): 15nmxl0nm. Tunnel current; 
40 pA. 
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FIG. 4: Color online. Schematic representatioiP^ of some prototypical geometries of OCHs-styrene (a,c,d) and CF3- styrene 
(b) molecules on H terminated Si(lOO). Si,C,0,F and H atoms are blue, black, red, green and white respectively. In each case 
the molecule is shown together with the two Si atoms of the Si dimer to which the molecule bonds and two Si atoms of an 
adjacent Si dimer row. In the T-tethered geometries (a,b and d) the molecule is located mainly over the trench between Si 
dimer rows while in the D-tethered geometry (c) the molecule is mainly over the Si dimer row to which it is bound. The OC 
bonds of the OCH3 groups are T-oriented (towards the trench) in (a) and (c) or D-oriented (towards the dimer to which the 
molecule binds) in (d). 
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FIG. 5: Color online. Calculated current / flowing between the tungsten STM tip and a CFs-styrene/OCHa-styrene molecular 
chain on silicon at some positive substrate biases vs. STM tip position along the chain at constant tip height. A side view of 
a part of the molecular chain around the CFs-styrene/OCHa-styrene junction is shown at the top of the figure. Atoms are 
colored as in Fig|4] A top view of the positions of some key atoms is shown in Fig. [sja) together with the trajectory of the 
tungsten STM tip for the calculated current profiles. Curve L (black) is for a very low bias for which the STM tip Fermi level 
is between the bottom of the silicon conduction band at the Si surface and the lowest energy state derived from the molecular 
LUMOs. Curve E (red) is for a somewhat higher (but still low) bias for which the STM tip Fermi level is just above the lowest 
energy state derived from the molecular LUMOs. Curve H (blue) is for a still higher bias for which the STM tip Fermi level 
is near the top of the band of states derived from the molecular LUMO. Black bullets (red diamonds) at the top indicate the 
lateral locations of the C atoms of the molecular CF3 groups (O atoms of the OCH3 groups). 
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FIG. 6: Color online. Calculated current / flowing between the tungsten STM tip and a CFa-styrene/OCHa-styrene molecular 
chain on silicon at some negative substrate biases (flUed state imaging) vs. STM tip position along the chain at constant tip 
height. The geometry of the chain is that in Fig. |5]and Fig. [sja) where the trajectory of the STM tip is also shown. Curve L 
(black) is for a very low bias for which the STM tip Fermi level is between the top of the silicon valence band at the Si surface 
and the highest energy state derived from the molecular HOMOs. Curve E (red) is for a somewhat higher but still low bias 
for which the STM tip Fermi level is just below the highest energy state derived from the molecular HOMOs. Curve H (blue) 
is for a still higher bias for which the STM tip Fermi level is just below the bottom of the band of states derived from the 
OCHs-styrene molecular HOMO. The black bullets (red diamonds) at the top of the plot indicate the lateral locations of the 
carbon atoms of the molecular CF3 groups (O atoms of the OCH3 groups). 
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FIG. 7: Calculated electrostatic shifts of atomic orbital energies E„ defined by Eqn. ([T]) vs. lateral position of the orbital 
along the molecular chain for various molecular chain geometries. Red curves (Ph) show the average of En over the six carbon 
atoms of the benzene ring of each molecule (scale on the left axis), black curves labelled Si (right axis) show E„ for the Si 
atoms to which the molecules bond. Blue vertical dashed lines separate the CFs-styrene part of the chain on the left from 
the OCH3- styrene on the right. Case (a) is for the molecular chain geometry discussed in Section IV the calculated current 
profiles for this geometry are shown in Figs |5] and |6] The horizontal black dotted line in (a) indicates En for surface Si atoms 
if the molecules are replaced by H atoms. The distinguishing structural features of the molecular chains with the electrostatic 
profiles shown in panels (a)-(e) are presented in Fig. |8|a)-(e), respectively. 



26 



%.«.«l.«,«.J» •••Jl ^ A % ••••••• w» ^ 

(a) ^ ^ ^ "5 ^ t ^ c 6 4i 4! t *r *r ^ 



(b) 
(c) 
(d) 



•f* f*<*^* A ••L*** ••••••• 

• At^.f^t***^*** •••••• 

•«.%«^«.%J*«%A^*.t % AAj* • 



FIG. 8: Atomic positions projected onto the Si (001) surface for some CFs-styrene/OCHs-styrene molecular chains on silicon 
studied theoretically in this work. Each panel shows two rows of Si atoms belonging to different Si dimer rows. The molecules 
are tethered to the lower row of Si atoms in each panel. The C atoms that bond to the Si atoms are orange. The C atoms that 
bond to the C atoms that bond to the Si atoms are violet. The C atoms of the CF3 and OCH3 groups are black. The O atoms 
are red and the F atoms are green. The C atoms belonging to the benzene rings that bond to O atoms or to C atoms not 
belonging to the benzene rings are white. The other other C atoms of the benzene rings, the H atoms and the other Si atoms 
are not shown for clarity. The dashed lines show the paths of the apex atom of the STM tip for the theoretical STM current 
profiles presented in this paper, (a) is the structure for Fig. [sjandjo] (b),(c),(d),(e) correspond to Figs 9|10|11|12 respectively. 



27 




-4-3-2-101234 



Position (nm) 



FIG. 9: Calculated current I flowing between the tungsten STM tip and a CFs-styrene/OCHa-styrene molecular chain on 
silicon vs. STM tip position along the chain at constant tip height for the molecular geometry shown at the top of the figure 
and in Fig. [8|b) where the trajectory of the STM tip is also shown. Notation as in Figjs] Curve L (black) is for a very low 
negative substrate bias for which the STM tip Fermi level is between the top of the silicon valence band at the Si surface and 
the highest energy state derived from the molecular HOMOs. Curve E (red) is for a somewhat higher but still low negative 
substrate bias for which the STM tip Fermi level is just below the highest energy state derived from the molecular HOMOs. 
Curve E' (violet) is for a low positive substrate bias for which the STM tip Fermi level is just above the lowest state derived 
from the molecular LUMO. 
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FIG. 10: Calculated current / flowing between the tungsten STM tip and a CFs-styrene/OCHs-styrene molecular chain on 
silicon vs. STM tip position along the chain at constant tip height for the molecular geometry shown at the top of the flgure 
and in Fig. |8|c) where the trajectory of the STM tip is also shown. The magenta lines mark the locations of the CF3- 
styrene/OCHs-styrene interface (left) and of the dislocation in the OCHs-styrene chain (right) in the image at the top of the 
flgure and in the current plots below. Notation as in Fig. 5: Curve L (black) is for a very low negative substrate bias for which 
the STM tip Fermi level is between the top of the silicon valence band at the Si surface and the highest energy state derived 
from the molecular HOMOs. Curve E (red) is for a higher but still low negative substrate bias bias for which the STM tip Fermi 
level is just below the highest energy state derived from the molecular HOMOs. Curve E' (violet) is for a positive substrate 
bias for which the STM tip Fermi level is just above the lowest state derived from the molecular LUMOs. 
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FIG. 11: Calculated current I flowing between the tungsten STM tip and a CFa-styrene/OCHa-styrene molecular chain on 
silicon vs. STM tip position along the chain at constant tip height for the molecular geometry shown in Fig. [8|d) where the 
trajectory of the STM tip is also shown. Curve E (red) is for a low negative substrate bias bias for which the STM tip Fermi 
level is just below the highest energy state derived from the molecular HOMOs. Curve E' (violet) is for a positive substrate 
bias for which the STM tip Fermi level is just above the lowest state derived from the molecular LUMOs. 
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FIG. 12: Calculated current I flowing between the tungsten STM tip and a CFs-styrene/OCHs-styrene molecular chain on 
silicon vs. STM tip position along the chain at constant tip height for the molecular geometry shown in Fig. [sje) where the 
trajectory of the STM tip is also shown. Curve E (red) is for a low negative substrate bias for which the STM tip Fermi level 
is just below the highest energy state derived from the molecular HOMOs. Curve E' (violet) is for a positive substrate bias for 
which the STM tip Fermi level is just above the lowest state derived from the molecular LUMOs. 
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FIG. 13: Color online, (a) Schematic top view (not to scale) of model single-triple CFs-styrene structure. Molecules are green, 
surface Si atoms blue, H atoms white. The single CFa-styrene row consists of 8 molecules, the triple row of 24 molecules. All 
molecules are T-tethered to Si; see Figgb). (b) Plots V, L and H are calculated current profiles for an STM tip trajectory 
at constant height along center row of CFa-styrene molecules in (a) for negative substrate bias. Plot V: Low bias; ti p F ermi 
level near highest Si valence band states. Plot L: Stronger bias but tip Fermi level above CFs-styrene HOMO energies.'^ Plot 
H: Still stronger bias; tip Fermi level well within the HOMO energy band of the single CFs-styrene row. As in experiment, 
contrast between triple and single CFs-styrene rows in blue profile L is much weaker than in red profile V. Dashed black curve: 
Calculated current profile along short CFa-styrene row on left in (a) at same tip height and bias as plot V. Black bullets locate 
C atoms of CF3 groups in the long CFs-styrene row. (c) Black solid curve (Si): Electrostatic electronic energy shifts En at Si 
atoms to which molecules of the long CFa-styrene row bond. Triple (single) rows are right (left) of vertical dotted line. Red 
curve (Ph): Average of _E„ over the six carbon atoms of the benzene ring of each molecule of the long CFs-styrene row (scale 
on right axis). Black dashed curve shows for comparison _B„ at Si atoms to which molecules of a single-row heterostructure of 
10 CFa-styrene molecules (on the left) and 10 OCHa-styrene molecules (on the right) bond, computed for a silicon substrate 
cluster having the same cross-section. 



